1288C - Two Arrays - CodeForces Solution

#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#include <functional>
using namespace std;
using namespace __gnu_pbds;
typedef tree<int, null_type, less<int>, rb_tree_tag,
             tree_order_statistics_node_update>
    ordered_set;
typedef long long int lli;
#define rep(i, a, n) for (lli i = a; i < n; i++)
#define rrep(i, n, a) for (lli i = n; i >= a; i--)
#define vlli vector<lli>
#define mii map<lli, lli>
#define umii unordered_map<lli, lli>
#define plii pair<lli, lli>
lli mod = 1e9 + 7;
void fastio()
{
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
}
lli power(lli a, lli b)
{
    lli res = 1;
    while (b)
    {
        if (b % 2 == 0)
        {
            a = (a % mod * a % mod) % mod;
            b = b / 2;
        }
        else
        {
            res = (res % mod * a % mod) % mod;
            b--;
        }
    }
    return res % mod;
}
lli mod_inv(lli a)
{
    return power(a,mod-2);
}
int main()
{
    fastio();
    int n,m;
    cin>>n>>m;
    lli a=1,b=1,c=1;
    rep(i,1,n+2*m)
    {
        a=(a*i)%mod;
        if(i<n){
            c=(c*i)%mod;
        }
        if(i<=2*m){
            b=(b*i)%mod;
        }
    }
    // cout<<a<<" "<<b<<" "<<c<<"\n";
    lli ans=a;
    ans=(ans*mod_inv(b))%mod;
    ans=(ans*mod_inv(c))%mod;
    cout<<ans;
    return 0;
}